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SUMMARY 


Three methods have been developed for estimating the damping in 
roll through the normal flight range of lift coefficient for wing plan 
forms having various sweep angles, aspect ratios, and taper ratios. 

Values of the damping in roll calculated by the three methods were 
compared with experimental values. The most complete method in which 
all known factors affecting the dancing in roll were considered 
appeared to give the best quantitative agreement with experiment for 
approximately 60 percent of the cases investigated. Another method, 
in which the value of the damping in roll at zero lift coefficient is 
modified only in accordance with variations in finite-span lift-curve 
slope, is a lmo st as reliable as the most complete method. The most 
important factor considered in the analysis therefore appears to be 
the variation of the finite-span lift-curve slope. 


INTRODUCTION 


Results of tests conducted on swept wings (reference 1, for 
example) have indicated that, although the present theories used to 
calculate the damping in roll (references 2 to 6) are quite accurate 
at the low lift coefficients, large discrepancies may exist at 
moderate and high lift coefficients. 

The present analysis was made to determine the factors that 
contribute to the discrepancies at high lift coefficients and to 
devise procedures by which improved estimates could be made. The 
procedures developed from the present analysis are based on simplified 
theoretical concepts and utilize measured lift and drag data. 
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Comparisons are presented between experimental values of the 
damping in roll and values estimated by three methods. 


SYMBOLS 


The symbols used in the analysis and in the presentation of the 
results are defined herein. 


°L 


lift coefficient 


Lift 

ipV2sj 

,2 


C D 


drag coefficient ( Dr a g \ 


L 


rolling-moment coefficient [ Rolling moment 

V ^pV 2 Sb 


'n 


yawing— moment coefficient ( Yawing moment 

JpV 2 Sb 


P 

S 

V 

b 

y 

a 

A 


mass density of air 
wing area 
flight velocity 

wing span, measured perpendicular to plane of symmetry 

spanwise distance from plane of symmetry to any section on 
wing quarter-chord line 

angle of attack, measured in plane of symmetry 

sweep angle of wing quarter-chord line, positive for 
sweepback 


aspect ratio (b 2 /!3) 
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pb/2V wing— tip helix angle 

p angular velocity about X-axis 

a Q section lift— curve slope 

c 2 section primary force coefficient (see reference 4) 

C = 

^ 8a 



Kl, K 2 constants that are functions of wing plan form 

0 = \fl — M^cos^A 

M 

Subscripts: 

1 induced 

o profile 

C L at any lift coefficient unless specified differently 


Mach number (Velocity o f free stream 
\ Velocity of sound , 


R 


right wing panel 
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L left wing panel 

M at any Mach number unless specified differently 


ANALYSIS 

Development of Methods 


The fact that the damping in roll is a function of both the 
lift-curve slope and the drag has been demonstrated by Glauert (refer- 
ence T). In most instances, for straight wings of moderately high 
aspect ratio, the lift— curve slope changes little through the lift- 
coefficient range and the drag contribution is relatively unimportant. 
Potential— flow values of Cj at zero lift coefficient, therefore, 

have generally been considered satisfactory at all lift coefficients 
below the stall. 


A brief summary of the methods for calculating 
presented in' references 2 to 6 is given herein. 


C 7 at zero lift . 
*P 


Charts presenting values of the damping in roll based on lifting- 
line theory at zero lift coefficient for unswept wings having various 
aspect ratios and taper ratios sire given in reference 2. These values 
were calculated from spanwise lift distributions that • correspond to the 
rolling motion. The results given, in reference 2 were corrected for 
lifting-6urface effects by application of an effective edge-velocity 
correction in reference 3* -An approximate method of modifying the 
results of reference 3 for the effects of sweep is presented in refer- 
ence 4. A similar, but more refined, method which accounts for the 
effects of sweep on the edge-velocity correction factor is given in 
reference 5. A more rigorous method of calculating C j for wings of 


arbitrary plan form at zero lift coefficient is presented in reference 6 
This method consists of an application of the theory of Weissinger 
(reference 8) for determining the additional span loading during roll. 


The values for the damping in roll presented in reference 6 are 
considered to be the most reliable and are used herein at zero lift 
coefficient.' The results from reference 6 have been extended to an 
aspect ratio- of 10 and are presented in figure 1. The values 


(Cl 


c T =o 


presented in figure 1 are for a section lift-curve slope 


of 2«. The methods of reference 4 provide a convenient basis, for 
correcting these results to any value of section lifi>-curve slope. 


As 
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indicated in reference 6, the correction is 



A ± 4 cos A 


2a 

a 


A + 4 cos A 


(1) 


For low-aspect— ratio or swept wings of the type considered for 
high-speed flight, the finite— span lift— curve slope may vary consider- 
ably through the lift-coefficient range and the drag at high lift 
coefficients may become large. A procedure based on potential— flow 
considerations for calculating Cj at high lift coefficients can 

hardly be expected, therefore, since such behavior is associated with 
separation of flow that results in important local changes in the 
characteristics. 

The method of reference 4 can be used to demonstrate that, to a 
first approximation, variations in the lift-curve slope will 
affect Cj p in the same proportion as . Cj^. Therefore, if the effects 

of drag are neglected, Ci can be written as 



\ ^ t 

An analysis of the drag contribution for straight wings of 
elliptic-chord distribution (reference 9) indicates that the increment 

of Cj . due to drag is equal to £ Cp. If the same increment can be 

applied to swept wings, the resulting expression for C is 

l P 



(3) 
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Ths problem of accounting for the drag of swept wings is . actually 
somewhat more complicated than indicated by equation (3) since the 
profile drag and induced drag are not of equal importance and, there- 
fore, must be considered separately. That is. 



in which the assumption is made that 


ClT 


and and K 2 are constants that are functions of the wing plan 

form. 

According to the strip— theory procedure of reference 4, the 
induced-drag contribution to Cj p can be derived as follows. 


The rolling moment due to induced drag for a constant-chord swept 
wing in roll can be expressed as 



■'b/2 


C2 t sin A tan A 

Xj 



c^ sin A tan A 



c 2t cos A ^ - C2 cos A 


y <*y 


(5) 
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where, as Indicated In reference 4, 


c 2 t 


1 

jtA cos A 



Aa Q cos A y 
A + 4 cos A b/2 



°2 R 



\ o 

E*'| 
27, 




Substituting the values c^ and C 2 R In equation performing 

the integration, differentiating with respect to pb/2V, and 
assuming = 2 n result in the following expression: 


(° z p) 


°L" 


A 


°Di 


f i + p s1h2a A + 2 cos A \ 

5 nA cos* A y 8ln A A + 4 cos AJ 


For zero sweep this expression reduces to 



(6) 


which is in agreement with the result given in reference 9 for constant- 
chord wings. It can be demonstrated that for swept wings having 
elliptic-chord distributions, the constant l/6 in equation (6) is 
replaced by 1/8, as was used in connection with the total-drag term 
in equation (3). The constant corresponding to the elliptic-chord 
distribution is regarded as being more reliable for wings of practical 
design, and therefore the induced-drag contribution to C 7 can be 

tp 
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written as 


(%), 


= - i ° L (l + 2 81I& A ± 2 ° 0S A 


°I>i 


8 nA cos^ A 


A + 4 cos A 


) 


(7) 


Values obtained from equation (7) are presented in figure 2 as a plot 


of 


( Cz p) 


°Di, 


Cjf against aspect ratio for various sweep angles. 


As in the case of the induced drag, the theoretical constant for 
the profile-drag contribution to for an elliptic-chord distribu- 

tion is l/8, which is the same as for straight wings. 

An equation for the damping in roll of swept wings, which accounts 
for the induced drag, the profile drag, and for variations in finite- 
span lift-curve slope, can therefore be written as 


C, = 
ip 


* \ , / c \ 
( C 'pVo(^*( Ci p) 


ODj 


— ^ Ct) 

8 Vc 


(8) 


An attempt was made to evaluate empirical values of Kg (equa- 
tion (4) ) by procedures similar to that used for C n in reference 1, 

P 

but the value of l/8 appeared to be as good as any that could be obtained 


Application of Methods 

From the foregoing analysis, three methods of calculating are 

indicated by equations (2), ( 3 ), and (8): 

Method 1 


/ \ Wqt 

p ' ( !p k=° (Svo 
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Method 2 




The values of 



used in methods 1, 2, and 3 are presented 


in figure 1 for zero Mach number and a section lift-curve slope of 2n. 
As was previously pointed out, the values from figure 1 can he corrected 
to any section lift-curve slope by means of the relation 



A + 4 cos A 



A + 4 cos A 


Methods of applying corrections for the effects of compressibility 
within the subcritical range of Mach number, are indicated in refer- ’ 
ences 6 and 10,. The method of reference 10, though less rigorous than 
that of reference 6, has been found to be about as reliable as and 
somewhat more convenient than the method of reference 6. According to 
reference 10, the correction for compressibility is 


(C, 

'p)c L =o 

_ A + 4 cos A 
Aj3 + 4 cos A 

( C!p )cl=c 


_ 

M 
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The ratio 



used in methods 1, 2, and 3 is 


defined as the ratio of the finite-span lift-curve slope at a parti- 
cular lift coefficient to the finite-span lift-curve slope at zero lift 
coefficient. This ratio is obtained from the measured lift character- 
istics of the particular plan form. In order to estimate values 
of Cj at high Mach numbers the lift data corresponding to the part- 


icular Mach number should be used. Similarly, in method 2, the drag 
data corresponding to the particular Mach number should be used in order 
to account for compressibility effects. 


The values of /C7 \ 

l *7% 


used in method 3 may be obtained from 


figure 2. Corrections to this increment for variations in section 
lifi>-curve slope or Mach number should not be necessary, since this 
increment is small except at high lift coefficients where the Mach 
number generally is low. 


The profile- drag term C Dq 


of method 3 is defined as 


n _ c _2Lf 
VDo”-T> *A 


where Cp and should be taken from experimental data at the 

particular Mach number in question. 


RESULTS AND DISCUSSION 


Experimental data of 



obtained for the wings listed in 


table I (taken from references 1, 11, and NACA tests) are compared in fig- 
ures 3 to 17 with the results calculated by each of the three methods. 

The experimental data of C L , Cp, and Cp used in the calculations 


are also presented. 


Examination of figures 3 to 17 shows, that method 1, in which the 
value of Cj at = 0 is modified only in accordance with 
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variations in Cj^, is almost as reliable as method. 3. The variation 

of finite-span lift-curve slope therefore appears to be the most 
important factor considered, in the analysis. 

Method 3 shows a considerable error for the swept wings of aspect 
ratio 1.34 (wings 6 and 7) presented in figures 5 and 6 . The error 
seems to result from the fact that the large values of fCi \ in— 

V lp JpDi 

dicated in figure 2 apparently are not realized for these wings. This 
result might be expected, since the assumptions regarding the distri- 
bution of forces on the wings could hardly be reliable at so low an 
aspect- ratio. 

Of these three methods, method 3, which includes separate 
considerations of the effects of induced and profile drag on C 1 , 

tp 

appears to give the best quantitative agreement with experiment for 
approximately 60 percent of the cases investigated. 

In the case of the two sveptforward wings presented in figure 15, 
the agreement between the experimental and calculated results is poor 
at the moderate and high lift coefficients. The large increase 
in C 1 with lift coefficient for these wings is not accompanied by 

large changes in the variation of Cj^ with lift coefficient. This 

condition probably accounts for the poor agreement since, as has been 
pointed out, the variation in the finite— span lift-curve slope is the 
most important factor considered in the analysis. 


CONCLUDING REMARKS 


Three methods have been developed for estimating the dampin g in 
roll through the normal flight range of lift coefficient for wing 
plan forms having various sweep angles, aspect ratios, and taper 
ratios. 

Values of the damping in roll calculated by the three methods 
were compared with experimental values. The most complete method 
(method 3) in which all known factors affecting the dampin g in roll 
were considered appeared to give the best quantitative agreement with 
experiment for approximately 60 percent of the cases investigated. 
Another method (method 1 ), in which the value of the damp ing in roll 
at zero lift coefficient is modified only in accordance with variations 
in finite-span lift-curve slope, is almost as reliable as the most 
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complete method. The most important factor considered in the analysis 
therefore appears to he the variation of the finite-6pan lift-curve 
slope. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va. , June 7 , 194-9 
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TABLE I 

SUMMARY OF PERTINENT INFORMATION REGARDING WINGS CONSIDERED IN ANALYSIS 


Wing 

Sweep angle, 
A 

(teg) 

Aspect ratio, 
A 

Taper ratio, 
X 

Test Reynolds 
number 

Airfoil section 

Figure 

Source 

1 

45 

2.61 

0.25 

1.56 X 10 6 

NACA 0012 

3 



2 

'45 

2.61 

.50 

1.45 

NACA 0012 

3 











> NACA tests 

3 

45 

2.61 

1.00 

1.39 

12 percent biconvex 

4 



4 

45 

2.61 

1.00 

1.39 

NACA 65^-012 

4 



5 

0 

1.34 

1.00 

1.99 

NACA 0012 

5 



6 

45 

1.34 

1.00 

1.97 

NACA 0012 

5 



7 

60 

1.34 

1.00 

1.97 

NACA 0012 

6 



8 

0 

2.61 

1.00 

1.39 

NACA 0012 

6 



9 

45 

2.61 

1.00 

1.39 

NACA 0012 

7 


\ Reference 1 

10 

60 

2.61 

1.00 

1.37 

NACA 0012 

7 



11 

0 

5.16 

1.00 

.98 

NACA 0012 

8 



12 

45 

5.16 

1.00 

• 97 

NACA 0012 

8 

/■ 


13 

60 

5.16 

1.00 

.76 

NACA 0012 

9 

✓ 


14 

2 

10.00 

.50 

.163 

Rhode St. Genese 33 

9 



15 

42 

5.90 

.50 

.231 

Rhode St. Genese 33 

10 


> Reference 11 

16 

62 

2.50 

.50 

.326 

Rhode St. Genese 33 

10 

j 


17 

0 

4.62 

.55 

5.6 

°NACA 0015 

11 



18 

30 

4.84 

.44 • 

6.20 

a NACA 0015 

11 



19 

45 

3.64 

.42 

7.80 

SflACA 0015 

12 



20 

3.6 

4.00 

.60 

.710 

NACA 65 AOO 6 

12 



21 

32.6 

4.00 

.60 

.710 

NACA 65 AOO 6 

13 


i 

22 

46.7 

4-. 00 • 

.60 

.710 

NACA 65 AOO 6 

13 



23 

36.9 

4.00 

0 

1.23 

NACA 0012 

14 










/ NACA tests 1 

24 

52.2 

2.31 

0 

1.62 

NACA 0012 

14 



25 

-45 

2.61 

1.00 

1.39 

NACA 0012 

15 • 



26 

^5 

3.12 

.38 

| 8.95 

a NACA 0015 

15 



27 

45 

4.00 

'.60 

I .710 

NACA 65 AOO 8 

16 



b 28 

45 

4.00 

.60 

.710 

NACA 65 AOO 8 

16 



29 

37.5 

3.00 

.49 

1.02 

NACA 23012 

17 



b 30 

37-5 

3.00 

.49 

1.02 

NACA 23012 

; x 7 

j 



Varied from NACA 0015 section at root to NACA 23 OO 9 section at tip. 
^Tested with full-span slats. 













Experimental 



Figure 3.- Estimation of the damping in roll for wings 1 and 2. 











Experimental 



Figure 6.— Estimation of the damping in roll for wings 7 and 8. 










Method 1 



Figure 9«— Estimation of the damping in roll for wings 13 and 14 










Figure 12.— Estimation of the damping in roll for wings 19 and 20. 







Method 2 



Figure 14.— Estimation of the damping in roll for vings 23 and 2k 


Method 1 



Figure 15 •— Estimation of the damping in roll for wings 25 and 26 


Experimental 



Figure 1 6 .— Estimation of the damping in roll for wings 27 and 28 





Figure 17 •— Estimation of the damping in roll for wings 29 and 30 




